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ABSTRACT 

We investigate whether the parsec-scale jets of quasars 3C 345 and 3C 279 are dom- 
inated by a normal (proton-electron) plasma or a pair (electron-positron) plasma. 
We first present a new method to compute the kinetic luminosity of a conical jet by 
using the core size observed at a single very long baseline interferometry frequency. 
The deduced kinetic luminosity gives electron densities of individual radio-emitting 
components as a function of the composition. We next constrain the electron density 
independently by using the theory of synchrotron self-absorption. Comparing the two 
densities, we can discriminate the composition. We then apply this procedure to the 
five components in the 3C 345 jet and find that they are pair-plasma dominated at 
14 epochs out of the total 19 epochs at which the turnover frequencies are reported, 
provided that the bulk Lorentz factor is less than 15 throughout the jet. We also in- 
vestigate the composition of the 3C 279 jet and demonstrate that its two components 
are likely pair-plasma dominated at all the four epochs, provided that their Doppler 
factors are less than 10, which are consistent with observations. The conclusions do 
not depend on the lower cutoff energy of radiating particles. 

Key words: galaxies: active — quasars: individual (3C 345) — quasars: individual 
(3C 279) — radio continuum: galaxies 
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1 INTRODUCTION 

The study of extragalactic jets on parsec scales is astrophysi- 
cally interesting in the context of the activities of the central 
engines of active galactic nuclei (AGN). In particular, a de- 
termination of their matter content would be an important 
step in the study of jet formation, propagation and emission. 
There are two main candidates for their matter content: a 
'normal plasma' consisting of protons and relativistic elec- 
trons (for numerical simulations of shock fronts in a very 
long baseline interferometry jet, see e.g., Gomez et al. 1993, 
1994a,b), and a 'pair plasma' consisting only of relativistic 
electrons and positrons (for theoretical studies of two-fluid 
concept, see Sol, Pelletier & Asseo 1989; Despringre & Fraix- 
Burnet 1997). Distinguishing between these possibilities is 
crucial for understanding the physical processes occurring 
close to the central engine (presumably a supermassive black 
hole) in the nucleus. 

Very long baseline interferometry (VLBI) is uniquely 
suited to the study of the matter content of pc-scale jets, 
because other observational techniques cannot image at mil- 
liarcsecond resolution and must resort to indirect means of 
studying the active nucleus. Recently, Reynolds et al. (1996) 



analyzed historical VLBI data of the M87 jet at 5 GHz 
(Pauliny-Toth et al. 1981) and concluded that the core is 
probably dominated by an e plasma. In the analysis, they 
utilized the standard theory of synchrotron self absorption 
(SSA) to constrain the magnetic field, B [G], and the proper 
number density of electrons, A'^^ [1/cm''], and derived the fol- 
lowing condition for the core to be optically thick for SSA: 
N^B^ > 25^ax, where ^max refers to the upper limit of the 
Doppler factor of the fluid's bulk motion. This condition is, 
however, applicable only for the VLBI observations of M87 
core at epochs September 1972 and March 1973. Therefore, 
in order to apply the analogous method to other AGN jets 
or to the M87 jet at other epochs, we must derive a more 
general condition. 



On these grounds, Hirotani et al. (1999, hereafter Paper 
I) generalized the condition N^B^ > 25,^ax and applied it to 
the 3C 279 jet on parsec scales. In that paper, they revealed 
that the core and components C3 and C4, of which spectra 
are reported, are dominated by a pair plasma. It is interest- 
ing to note that the same conclusion that component CW in 
the 3C 279 jet is dominated by a pair plasma was derived by 
an independent method by Wardle et al. (1998), who stud- 
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ied the circularly and linearly polarized radio emission from 
3C 279. 

Subsequently, Hirotani et al. (2000, hereafter Paper 
II) applied the same method to the 3C 345 jet. Deducing 
the kinetic luminosity from a reported core-position offset 
(Lobanov 1998), they demonstrated that components C2, 
C3, and C4 at epoch 1982.0, C5 at 1990.55, and C7 at four 
epochs are dominated by a pair plasma. 

In the present paper, we revisit the two blazers, adding 
spectral information at 2 epochs for 3C 279 and 11 epochs 
for 3C 345. In the next section, we first improve the method 
presented in Paper II and demonstrate that the kinetic lu- 
minosity, Z/kin, can be inferred from the VLBI core size at a 
single frequency, rather than from the core-position offset, 
which requires more than two VLBI frequencies. We then 
provide the method to deduce the electron density, A^^* (nml), 
from I/kin, assuming a normal-plasma dominance. In §3, we 
describe an independent method to deduce the electron den- 
sity, A'^(SSA), by utilizing the theory of SSA; this method 
is based on the generalized condition of N* B^ > 2(5max~^ in 
Reynolds et al. (1996). If iVe (nml) < Af*(SSA) holds for a 
radio-emitting component, we can rule out the possibility of 
a normal-plasma dominance. In § 4, we apply the method 
to the 3C 345 jet and find that components C2, C3, C4, 
C5, and C7 are dominated by a pair plasma as a whole. We 
next investigate the 3C 279 jet in § 5 and demonstrate that 
components C3 and C4 are pair-dominated, provided that 
their Doppler factors are less than 10. In the final section, 
we discuss the possibility of the entrainment of the ambient 
matter as a jet propagates downstream. 

We use a Hubble constant Hq — 65h km/s/Mpc and 
go = 0.5 throughout this paper. Spectral index a is defined 
such that Sv oc v°' . 



2 KINETIC LUMINOSITY DEDUCED FROM 
THE CORE SIZE 

In this section, we present a new method to deduce the ki- 
netic luminosity of a conical jet from the core size observed 
at a single VLBI frequency in 



2.1-2.5. We then describe 



how to cor nput e the electron density from the kinetic lumi- 
nosity in § p.q . 

2.1 Conical Jet Geometry 

We assume that the parsec-scale jet close to the central en- 
gine propagates conically with a half opening angle % (fig-|l[)- 
Then the optical depth r for synchrotron self absorption at 
distance p from the injection point, is given by 



r„{p) = 



psmx 



+ ■ 



psmx 



sin(i^ + x) sin(i^ - x) 



(1) 



where ip is the viewing angle and a,, [1/cm] refers to the 
effective absorption (i.e., absorption minus stimulated emis- 
sion) coefficient. For a small half opening angle (x ^ 1), 
this equation can be approximated as 



Y 

r„{p) = 2p- a, 

simp 



(2) 



Note that x/ sin p denotes the projected, observed opening 
angle. 



injection 
piont 




Figure 1. Schematic figure of a conical jet in the core region 
with half opening angle x in the observer's frame. 



Since r^ and px are Lorentz invariants, we obtain 



^^ = ^i;^, (3) 

sm p sin p* 

where a quantity with an asterisk is measured in the co- 
moving frame, while that without an asterisk in the ob- 
server's frame. Since uai, is also Lorentz invariant, equation 
(|) gives 

!i^ = Ji = -^. (4) 

sin p V* 1 + z 

Combining equations (S) and (H), we obtain 



z 2px 



Dl 



(5) 



S sin(/5 " 6 smp{l + z)^ 

where Z)l/(1 + z)'^ is the angular diameter distance to the 
AGN, and 8d is the angular diameter of the component in 
the perpendicular direction of the jet propagation in radian 
units; the Doppler factor, S, is defined by 

(6) 



r(l-/3cOSV9A)' 



where F = l/i/l — /3^ is the bulk Lorentz factor of the jet 
component moving with velocity /3c, and pA refers to the 
viewing angle in the AGN-rest frame. It is related with p, 
viewing angle in the observer's frame, as 

cos (p — /3agn 



COS(y3A 

1 — PAGN COS (p 

where /3agn is defined by 



1 -I- /3agn 
1 — /3agn 



1 + z 



(7) 



(8) 



For 3C 345(z = 0.595) for instance, if a photon is ejected 
with (pA = 2° to the observer in the AGN-rest frame, we 
observe it propagates in the direction p — 3.2° due to the 
redshift (/3agn < 0). We should distinguish p and pi a, be- 
cause 5, which appears in equation (M), should be defined in 
the AGN-rest frame. 

We assume that the electron number density between 
energies ^rricC? and (7 -I- dj)mcC^ is expressed by a power 
law as 

dN* 

— ^ = fc*7^°'"^ (7,nin < 7 < 7max). (9) 
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Integrating dN*/d'y from 7min to 7max, and assuming 
7max ^ 7min and Q < 0, we obtain the electron number 
density 



n: = 



In 



-2a 






(10) 



Moreover, averaging over pitch angles of the isotropic elec- 
tron power-law distribution (eq. |9|), we can write down a^ 
in the co- moving frame of the jet as (Le Roux 1961, Ginzburg 
& Syrovatskii 1965) 



a* = C{a)ro k\ 






(_2a+3)/2 



(11) 



c/[eV(moc2)] = 1.063 X 10^'' GHz and 



where Uo = c/r, 
vb = eB/{2nmcc); e, me, and c refer to the charge on 
an electron, the rest mass of an electron, and the speed of 
light, respectively. The coefficient C{a) is given in Table 1 
of Gould (1979) and also in table 1 in this paper. 



2.2 Scaling Law 

Let us introduce a dimensionless variable r = p/ri to mea- 
sure p in parsec units, where ri — 1 pc. Assuming that there 
is neither particle production nor annihilation, we obtain the 
following scaling law for a conical jet from the conservation 
of the particles: 



iV* = Nir- 



(12) 



where n = 2; Ni refers to the value of A'^ at r = 1 (i.e., 
at 1 pc from the core in the observer's frame). In the same 
manner, the magnetic flux conservation in a stationary, su- 
personic magnetohydrodynamic jet leads to 



B = Bir 



(13) 



where m = 1; Bi refers to the values of B at r = 1. We 
further introduce the following dimensionless variables: 

2 , 



XB = l^Bi/l^o — 



riTo Ni, 
eBi 



27rmcC 



(14) 



Utilizing equation (llll) , and rewriting k* and z/b in terms of 
xn and xb, we obtain from the left equality in equation (BI) 



= C{a) 



2X 



sm <fi 7m 
a;Na;B 



^mi^) 



(15) 



where e = 3/2 — a. 

At a given frequency i/, the flux density will peak at 
the position where r^ becomes unity. Thus setting t — 1 
and solving equation (113) for r, we obtain the distance from 
the VLBI core observed at frequency v from the injection 
point (or the central engine) as (Lobanov 1998) 



r(u) = (xb''''F— 
V u 

where 

2x 



i/fcp 



/(«) 



fch 



C(a) 



sm if 7mi 



-2a ( 

. 2q \T 
:iiii \ r 



1 1/(^+1) 



Sn 



2a 



2a' 



(16) 

(17) 
(18) 
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(19) 



2.3 Core-Position Offset 

If we mesure rUA at two different frequencies (say v^, and 
Vh), equation (y_6|) gives the dimensionless, projected dis- 
tance of r{vz) — r{vh) as 



Arp 



(xB^-fUo)'^'"'^ 



l/fcr 



1/fcr 1/fer 



■ sm f. 



(20) 



Here, Arproj is in pc units; therefore, nArproj represents 
the projected distance of the two VLBI cores (in cm, say). 
That is, ri Arproj equals 4.85-10~^Armas£'L/(l-l-z)^ in equa- 
tion (4) in Lobanov (1998), where Armas refers to the core- 
position difference in mas. Defining the core-position offset 
as (Lobanov 1998) 



i/fe,- i/fc. 






;riz 



' P-^OJ l/fc. 



l/fcr 



we obtain 

^ = (4v-o 



,1/fcr 



smiy3 



(21) 



(22) 



To express xb in terms of a;N and ^rv , we solve equation 
(E2) for xb to obtain 



Xb = 



ri sin (f 



fcr/fcb 



(/^o 



-1/fcb 



(23) 



Note that xn is included in / = f(a). 

Setting fb — * oo in equation (|20[), we obtain the absolute 
distance of the VLBI core measured at f from the central 
engine as (Lobanov 1998) 



(^) 



-1/fcr 



ri sm if 



(24) 



That is, once fin/ is obtained from multi-frequency VLBI 
observations, we can deduce the distance of the synchrotron- 
self-absorbing VLBI core from the central engine, assuming 
the scaling laws of A'^ and B as equations ( |l2[ ) and ([l3|). 

We next represent xn and xb (or equivalently, Ni and 
Bi) as a function of fln^. To this end, we relate N* and B 
as follows: 



i\r* 2 



K 



B^ 

S-k' 



(25) 



When an energy equipartition between the radiating parti 
cles and the magnetic field holds, we obtain 

7niin 



K ^ 



il-)" 



(26) 



where the averaged electron Lorentz factor, (7_) becomes 



(7-) 



7 • fee 7 d'y 



fc*7^" ^d-y 



7mi 

2a 



\ O'ni ax / 'T^m i n J 



2a + l^'"'""(7max/7min)2--l ' 

for a < 0. In the special case a -^ —0.5, it reduces to 



(27) 
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(7-) = 7mm In 



O'max 



(28) 



Substituting N^ — Nir^^ and B — B\r^^ into equa- 
tion (ESI), and replacing A^i and Bi with a;N and x^, we 
obtain 

xn = xb (29) 

^ 7mm ^o 

It is noteworthy that the assumptions of n = 2 and m — 1, 
which results in fc = 1, in consistent with the energy 
equipartition at an arbitrary distance, r. Combining equa- 
tions (bs) and (E9|), we obtain 



Xb 



1 Lriy / l^Q 

ri sin ip 
X 



(5-2a)/(7-2a) 



7rC(a 



K n 



-2a 



sm ip 7miii r-o 7mii 



1 + z 



-2/(7-2a) 



(30) 



For ordinary values of a(< 0), sb decreases with increasing 
K, as expected. The particle number density, xn, can be 
readily computed from equation (Ea). 



2.4 Kinetic Luminosity 

When the jet has a perpendicular half width R±_ at a certain 
position, Lkin and A'^^ are related by 

Lkin = -kRa^^Pc-VN: ■ (r- 1) {{'y-)m,c^ + {7+)m+c^) , (31) 

where {7+) refers to the averaged Lorentz factors of posi- 
tively charged particles, and m+ designates the mass of the 
positive charge. 

For a conical geometry, we can put R±_ — px- Substi- 
tuting this relation into equation (pi|), and replacing N* p^ 
with 



N*p = iVin = — -a;N 
ro 

■R K r^^ 2 
- tXb , 



^ 'Yniin ^o 



we obtain 



Lkin = c\inK-^pr{r - 1) 



ro 



1 X^ru 



(32) 



2{5-2q)/(7-2q) 



7rC(a) K ri -2a 



sm p 7min ro 7mi 



2a / d \ 
^ [l + zj 



Vo r\ sm ip 

3/2-a-| -4/(7-2<.) 



(33) 



we have Ckin = 7i"^{7-)'7icC^/7min for a pair plasma, while 
Ckin = vr^mpC''/(27min) for a normal plasma, where rrip 
refers to the rest mass of a proton. It should be noted that 
7min for a pair plasma takes a different value from that for 
a normal plasma. 



2.5 Jet Opening Angle 

Let us now consider the half opening angle of the j et, and 
further reduce the expression of Lkin. As noted in 
obtain 



2.4, we 



Ra_ ^ 0.5ed,co.c(t^) I>L 

P ^ P (1 + - 



(34) 



for a conical geometry, where Sd,corc(i^) is the angular diam- 
eter of the VLBI core at a frequency v in the perpendicu- 
lar direction to the jet-propagation direction on a projected 
plane (i.e., a VLBI map). The luminosity distance is given 
by 



Dl = 4.61 X lO^fe- Vi ^°" + ('^'^ ' ^)("^ + ^^^^^^^+1. (35) 
On the other hand, equation (EJ) gives 



or\ 



-l/fcr 



svnp 



(36) 



Substituting equations (p6|), (pq) into 



we obtain 



1 x^rv _ igo£+ (go - 1)(-1 + \/2go2~fT) 



vq r\ sm p 



90^(1 + z) = 



i/fc. 



fo 



(37) 



As a result, the factor containing x and Sl,.j^ in equation 
(B3) can be simply expressed in terms of Sd,corc(j^)- That is, 
to evaluate Lkin, we do not have to stop on the way at the 
computation of VLrv, which requires more than two VLBI 
frequencies. We only need the VLBI-core size at a single 
frequency, 6'd,corc(i^), in addition to F, i/p, a, 7min, z, K, and 
Ckin- The information of the composition is included in Ckin- 
If we have to know x, we additionally need Q,rv 



2.6 Electron density deduced from kinetic 
luminosity 

Once Lkin is known by the method described above or by 
some other methods, we can compute the electron number 
density as a function of the composition. Solving equation 
(|l|) for iV*, and utilizing R^ = (6'd/2)DL/(l+z)^ we obtain 
for a normal plasma (i.e., for {7+) = 1 and m+ = rrip) 



iVc*(nml) = 



5.93 X 10~^/i^ 



r(r- 

go'(i 



-A 



mas/ 



zqo + (go - i)(-i -F V2go2 + 1) 



1 + (7-)me/mp 



(38) 



where L46 — Lkin/ (10*'' ergs s~^). If A'^(nml) becomes less 
than the electron density deduced independently from the 
theory of synchrotron self-absorption (SSA), the possibility 
of a normal plasma dominance can be ruled out. 

It is worth comparing Af*(nml) with A'^(pair), the par- 
ticle density in a pure pair plasma. If we compute Lkin by 



using equation 

N* (pair) _ 
iVe*(nml) 



033 



and (pTJ), we obtain the density ratio 



Lkin(pair) nipC^ + {7_(nml))mcC^ 



Lkin(nml) 
7min(nml) 
7min(pair) '' 



2(7- (pair) )mc 



(39) 



where 7min(pair) and 7niin(nml) refer to 7niin when a pair 
and a normal plasma dominates, respectively, (7- (pair)) and 
(7_(nml)) refer to the averaged Lorentz factor of electrons 
in a pair and a normal plasma, respectively. Therefore, if 
7min(pair) is comparable with 7niin(nml) ~ 100, there is lit- 
tle difference between A'^^(pair) and A'^c^(nml). 
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In short, if A'g*(nml) becomes much less than the A'^^ 
deduced independently by another method, we can rule out 
the possibility of a normal-plasma dominance and that of a 
pair-plasma dominance with 7min ~ 100 or greater. In the 
next section, we investigate how to deduce N* independently 
by the theory of SSA. 



3 SYNCHROTRON SELF-ABSORPTION 
CONSTRAINTS 

3.1 Magnetic Field Strength 

In this paper, we model a jet component as a homogeneous 
sphere of angular diameter 9d, containing a tangled mag- 
netic field B [G] and relativistic electrons which give a syn- 
chrotron spectrum with optically thin index a and maximum 
flux density Sm at the peak frequency v^. We can then com- 
pute the magnetic field density as follows (Appendix) : 



«-»-'"<')(<i:)"(^)'(t)^ 



(40) 



The coefficient b[a) is computed for a uniform sphere of 
plasma with angular diameter 9a,; their values are given in 
table 1. 



3.2 Proper Electron Density 

Substituting equation (hit) into (pi), and assuming 7min ^ 
7max, we obtain 



n:b- 



rriaC 






T-^(Q)7min^" 

C(q) -2a 



_ sin if {1 + z)^ / 1 + z 



m' 



-a + 2.5 



e, D. K s J • ■ (^^) 

Evaluating u at Um, and combining with equation (kol), we 
obtain (Marscher 1983) 

7min^'^ h(l + z)^qo^smip 



N:{SSA) = e(Q) 



-2a zqo + (go - 1)(-1 + V^qoz + 1) 



(I± 



KGUzJ 



^'^-s ^ S^ \ "'"+' / 5 - 2«-3 



Jy 



l + z 



(42) 



Vmas 

where 

e(a) = 1.71 X 10"'' x [2.79 x 10-*b(a)]"-''' x -^ (43) 

C(a) 

It looks like from equation (WJ) that A^'^ (SSA) depends 
on ^d and Vm very strongly. In the case of a = —0.75, for 
instance, we obtain A'c*(SSA) oc 5d~^°J^m~*. Nevertheless, 
as a radio-emitting components evolves along the jet, its 
dd increases due to expansion while its Um decreases due 
to (synchrotron-fadiabatic) cooling. As a result, these two 
effects partially cancle each other and reduce the variations 
of Af*(SSA) along the jet. 



3.3 Summary of the Method 

Let us summarize the main points that have been made in 
§§ 2 and 3. 

(1) Observing the core size, 6d,corc(i^), of a conical jet at a 
single frequency, v, we obtain 



z^o ri sin ip 



oc ^d.corc ■ I' 



l/fcr 



(44) 



by equation (Km . 

(2) Substituting equation (BTt) into (pq), we obtain the ki- 
netic luminosity, Lkin, in terms of Sd.corc ■ v^'^^ , a, F, ^, K 
and 7min for the core. Note that 7min takes difi^erent values 
between a normal plasma and a pair plasma. 

(3) Assuming a normal-plasma dominance, we obtain the 
proper electron density, A^^ (nml) from Lkin by equation ( p8| ) . 

(4) We can independently deduce A^* from the theory of 
SSA as AT* (SSA) by equation (||). 

(5) If A'^,*(nml) < A'^,*(SSA) holds, we can rule out the 
possibility of a normal-plasma dominance and that of a 
pair-plasma dominance with 7min ~ 100 or greater. That 
is, A'^^ (nml) <^ N^ (SSA) indicates the dominance of a pair 
plasma with 7min ^ 100. 



4 APPLICATION TO THE 3C 345 JET 

Let us apply the method described in the previous sec- 
tions to the radio-emitting components in the 3C 345 jet 
and investigate the composition. This quasar {z = 0.595; 
Hewitt & Burbidge 1993) is one of the best studied ob- 
jects showing structural and spectral variabilities on parsec 
scales around the compact unresolved core (for a review, 
see e.g., Zensus 1997). At this distance, 1 mas corresponds 
to 5.42/i~^ pc, and 1 mas yr~^ to /3app — 30.3^~^, where 
Ho = 65/ikms-^Mpc-\ 



4.1 Kinetic luminosity 

To deduce Lkin, we first consider 6d.corai^/i^o, where kr = 1 
is adopted in equation (pj). From the reported core size at 
22.2 GHz at 6 epochs by Zensus et al. (1995), at 5 epochs 
by Unwin et al. (1997), and at 3 epochs by Ros et al. (2000), 
we can deduce the averaged core size, Sd,coro, as 0.296 mas 
(fig. pi); here, we evaluate the core size with l.sVab for the 
latter three epochs, where a and b refer to the major and 
minor axes at the FWHM of the elliptical Gaussian pre- 
sented in Ros et al. (2000). From ^d,corc = 0.296 mas, we 
obtain Od,coici^/i^o = 6.18 x lO^^** as the averaged value 
over the 14 epochs at 22 GHz. In the same manner, we ob- 
tain the following averaged values at different frequencies: 
^d,corc = 0.258 mas and Qd,co-ccV jvf, = 3.64 x lO"'^* over 3 
epochs at 15 GHz (Ros et al. 2000); 6'd,corc = 0.436 mas and 
9d,corav/yo = 4.38 x 10"" over 7 epochs at 10.7 GHz (1 
epoch from Unwin et al. 1994; 4 epochs from Zensus et al. 
1995; 2 epochs from Gabuzda et al. 1999); 6'd,corc = 0.343 
mas and Od,corcv/i'o = 2.71 x 10~^* over 7 epochs at 8.4 GHz 
(1 epoch from Unwin et al. 1994; 3 epochs from Zensus et 
al. 1995; 3 epochs from Ros et al. 2000); &d,coro = 0.423 mas 
and Od,coTcv/vQ = 1.98 x 10"^* over 11 epochs at 5.0 GHz (4 
epochs from Brown et al. 1994; 1 epoch from Unwin et al. 
1994; 3 epochs from Zensus et al. 1995; 3 epochs from Ros 
et al. 2000). Taking the weighted average of 9d,coTcv/iyo over 
the 42 epochs at the different 5 frequencies, we obtain 

^dxcci^ =4.01x10-". 



1^0 



(45) 



We next consider the spectral index that reflects the 
energy distribution of the power-law electrons in the core. 
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Table 1. Table of constants 



a 





-0.25 


-0.50 


-0.75 


-1.00 


-1.25 


-1.50 


-1.75 


-2.00 


a 


0.2833 


0.149 


0.103 


0.0831 


0.0740 


0.0711 


0.0725 


0.0776 


0.0865 


C 


1.191 


1.23 


1.39 


1.67 


2.09 


2.72 


3.67 


5.09 


7.23 


rn>(0) 





0.187 


0.354 


0.503 


0.639 


0.762 


0.876 


0.981 


1.079 


(rm) 





0.165 


0.310 


0.440 


0.557 


0.663 


0.760 


0.848 


0.930 


6 





1.63 


2.33 


2.30 


2.05 


1.75 


1.48 


1.25 


1.05 



We thus obtain 



/ />/ . \ — 1.56 

Lkin = 4.68 X 10«ft-i" (^) K"-'' 



xr(r-i) 



ergs s (49) 



sin if \1 + z , 

for a normal plasma, that is, Ckin = 7r^mpC'^/(27niin) in 
equation (pd). 



Core size of 3C345 jet at 22 GHz 




1985 



Figure 2. Angular core size of 3C 345 at 22.2 GHz. The hori- 
zontal line represents the bet fit value. 



From infrared to optical observations (J,H,K,L bands from 
Sitko et al. 1982; J,H,K bands from Neugebauer et al. 1982; 
K band from Impey et al. 1982; IRAS 12, 25, 60, 100 microns 
from Moshir et al. 1990), we obtain 



Q = -1.15 



(46) 



for the core. This value is consistent with that adopted by 
Zensus et al. (1995). 

Using equations (Usf) and (M), we obtain from equa- 
tion (H) 



1 X^rv 



2(5-2a)/(7-2c«) 



fo ri sm (fi 
Moreover, a — —1.15 gives 



4.97 X 10"^^/i"^". (47) 



7rC(a) 



-2aK ri 



-| -4/(7-2a) 



1 oi ^, in-14r^-0.43 -0.5 

1.31 X 10 K 7mm 



(48) 



4.2 Kinematics near the core 

To constrain Lkin further, we must consider the kinematics 
of the core and infer F and tp. If the apparent component 
velocity, /3app, reflects the fluid velocity, the lower bound of 
F can be obtained by Fmin = iV/Sapp^ + 1- 

In the case of the 3C 345 jet, components are found 
to be 'accelerated' as they propagate away from the core. 
For example, for component C4, /9app/i < 8 held when the 
distance from the core was less than 1 mas (Zensus et al. 
1995); however, it looked to be "accelerated' from this dis- 
tance (or epoch 1986.0) and attained /3app^ ~ 45 at epoch 
1992.5 (Lobanov & Zensus 1994). It is, however, unlikely 
that the fluid bulk Lorentz factor increases so much during 
the propagation. If /3app reflected the fluid bulk motion of 
C4, the bulk Lorentz factor, F, had to be at least 45 until 
1992.5. Then /3app ~ 6 during 1981-1985 indicates that the 
viewing angle should be less than 0.05°, which is probably 
much less than the opening angle. Analogous 'acceleration' 
was reported also for C3 (Zensus et al. 1995) and C7 (Unwin 
et al. 1997). 

On these grounds, we consider the apparent 'accelera- 
tion' of the components in their later stage of propagation on 
VLBI scales are due to geometrical effects such as 'scissors 
effects' (Hardee & Norman 1989; Fraix-Burnet 1990), which 
might be created at the intersection of a pair of shock waves. 
Such shock waves may be produced, for example, by the 
interaction between relativistic pair-plasma beam and the 
ambient, non-relativistic normal-plasma wind (Pelletier & 
Roland 1989; Pelletier & Sol 1992). In their two-fluid model, 
the pair beam could be destructed at a certain distance from 
the core on VLBI scales due to the generation of Langmuir 
waves. 

Then, how should we constrain F and (fi for the 3C 345 
jet? In this paper, we assume that the radio-emitting com- 
ponents before their rapid acceleration represent (or mimic) 
the fluid bulk motion. Under this assumption, Steffen et 
al. (1995) fitted the motions of C4 (during 1980-1986) and 
C5 (during 1983-1989) when the components are within 2 
mas from the core by a helical model. They found F = 5.8 
for C4 and 4.6 for C5 with ip = 6.8°, assuming h = 1.54 
(i.e.. Ho ~ lOOkm/s/Mpc) . These Lorentz factors will be 
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~ 9 and ~ 7 for C4 and C5, respectively, if /i = 1 (i.e., 
Ho = 65km/s/Mpc). Subsequently, Qian et al. (1996) inves- 
tigated the intrinsic evolution of C4 under the kinematics of 
r = 5.6 and ip = 2° - 8°, assuming h= 1.54; the range of 
viewing angles is in good agreement with ip = 2° (for C4) 
- 4° (for C2) at 1982.0 obtained by Zensus et al. (1995), 
who used a larger value of F = 10 under h = 1.54, which 
corresponds to F ~ 15 for h = 1. 

In the present paper, we assume that the helical model 
describes the fluid bulk motion when the components are 
within 2 mas from the core and adopt F = 15 as the typical 
value. If F is less than 15, the possibility of a normal- plasma 
dominance further decreases. Close to the core, we adopt 
(/3 = 2°, because the viewing angle is suggested to decrease 
with decreasing distance from the core (Zensus et al. 1995; 
Unwin et al. 1997). If tfi is less than 2° in the core, Lkin, and 
hence A'^(nml) further decreases; that is, the normal-plasma 
dominance can be further ruled out. For F = 15 and ip = 2° 
(or equivalently tpA = 1.3°), we obtain 5 = 27.0 for the core 
and hence 



Lkin = 9.2 X l(f''h-^-^'^K' 



0.57 / 7min \ " 

V 100/ 



ergs s 



(50) 



assuming a normal plasma. Substituting equation (pOl) into 
(p8|), we obtain TV^ (nml) for each component having angular 
diameter 9d- 



4.3 Composition of individual components 

We can now investigate the composition of the radio- 
emitting components in the 3C 345 jet, using their spec- 
tral information. We assume that the jet is neither acceler- 
ated nor decelerated and apply a constant Lorentz factor of 
F = 15 for all the components, from the reasons described 
just above. To compute A'^* (SSA), we further need ip for each 
component at each epoch. (Note that ^ = 2° is assumed for 
the core, not for the components.) For this purpose, we as- 
sume that /3app reflects the fluid motion if /3app < VF^ — 1 
and that tp is flxed after /3app exceeds VF^ — 1. That is, we 
compute tp (or equivalently, i/3a) from 

^Papp 



tany^A 



/3app2 +5^-1 



if /?a 



Vr^ 



VF2-1 



if /3app > \/V^-l 



(51) 



The constancy of p (or i^^a) when /3app > VF^ — 1 — 15.0 
could be justified by the helical model of Steffen et al. (1995) , 
who revealed that the viewing angle does not vary so much 
after the component propagates a certain distance (about 1 
mas in the cases of components C4 and C5) from the core. 
The synchrotron self-absorption spectra of individual 
components are reported in several papers. They were first 
reported in Unwin et al. (1994), who presented (i^m,Sm,a) 
and 6d (or ^ in their notation) of components C4 and C5 
at epoch 1990.7. Subsequently, Zensus et al. (1995) gave 
{v^,Sm,ct), 6a, and /3app/i for C2, C3, and C4 at 1982.0. 
However, the errors are given only for 6^; therefore, we eval- 
uate the errors of the output parameters, 5, ip, B, N*{nra\), 
and A'^^(nml)/A^c (SSA), using the errors in 6d alone for 
C2, C3, and C4 at 1982.0. Later, Unwin et al. (1997) pre- 
sented (i'm,5'm,Q), ^d, and /9app/i of C5 (at 1990.55) and C7 
(at 1992.05, 1992.67, 1993.19, and 1993.55). More recently, 
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Lobanov and Zensus (1999) presented a comprehensive data 
of (!/m,5'm,Q) of C3, C4, and C5 at various epochs. We uti- 
lize Oa = 0.114 + 0.0658(epoch- 1979.50) for C4 (Qian et al. 
1996), Oa = 0.114 + 0.0645(epoch - 1980.00) for C5, which 
is obtained from the data given in Zensus et al. (1995) and 
Unwin et al. (1997), and Od = 0.09 -I- 0.45(psinv3/mas) for 
C3 (Biretta et al. 1986), where p sirup refers to the projected 
distance from the core. 

The results for individual components at various epochs 
are given in table 2 for K — 1. It follows that the ra- 
tio, A'^o*(nnil)/A'^^(SSA) becomes less than 1 for acceptable 
strength of the magnetic field, that is, ImG < B < lOOmG, 
say. For C4, we cannot rule out the possibility of a normal- 
plasma dominance at epochs 1985.8 and 1988.2. Neverthe- 
less, at these two epochs, the input parameters {vm, Sm, Q, 
9d) give too large B (> 380mG) for a moderate 5{< 30). For 
C5 at epoch 1984.2, the input parameters give too small B 
{— 9pG) for a moderate S{— 29.1). On these grounds, we 
discard these three cases. 

We present the ratio, A'^^(nml)/A'^^(SSA) as a func- 
tion of psimp in figure M. The solid line represents A'^^ — 
Ai'*(nml). It should be noted that Af!'(nml) is computed 
under {7-) = in equation (pSl). Thus, the solid line 
gives the conservative upper limit of Ai'*(nml)/Ai'*(SSA). In 
another word, if A'^(nml)/A'^(SSA) slightly exceeds unity 
(1.5, say), a normal-plasma dominance can be, in general, 
ruled out. There is also depicted a dotted line representing 
Nc = Ng{pa.iT). If the ratio appears under the dotted line, 
it indicates that Lkin is underestimated or that the input 
parameters give too small B and hence too large A'^(SSA). 

It follows from the figure that neither a normal-plasma 
dominance nor the pair-plasma dominance with 7min > 100 
are allowed within one-cr errors for the 14 epochs out of the 
19 epochs at which input parameters give ordinary B. The 
ratio shows no dependence on p under the assumption of 
constant F. 

It is 

noteworthy that N* {ami) /N^ (SSA) oc (7min/100)-^"-^'^'' 
for the spectral index of —1.15 for the core. Therefore, for 
a typical spectral index a ~ —0.75 for the jet components, 
the dependence on 7min virtually vanishes. That is, the con- 
clusion does not depend on the assumed value of 7min ~ 100 
for a normal plasma. 

If we assume instead F = 20 (rather than F — 15) and 
ip — 2° , we obtain the following kinetic luminosity 



1.3 X 10*" /i" 



-1.57 / 7min \ ' 

V 100/ 



K" 



ergs s 



(52) 



Evaluatin g th e viewing angles of individual components by 
equation ( pl| ) with F — 20, and excluding the three epoch 
data (C4 at 1985.8 and 1988.2, and C5 at 1984.2) that give 
extraordinary B, we can compute A'^^(nml)/A'^^(SSA) as pre- 
sented in flgure y. It follows from this figure that we cannot 
rule out the possibility of a normal-plasma dominance for 
such large Lorentz factors. The greater F we assume, the 
greater becomes A'"c*(nml)/A'"^(SSA). 

On these grounds, we can conclude that the jet is domi- 
nated by a pair plasma with 7min ^ 100 at 14 epochs out of 
the 19 epochs at which input parameters give extraordinary 
B, if we assume F < 15 throughout the jet. 
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Figure 3. The ratio Af*(nml)/Af*(SSA) as a function of the 
distance from the core. Above the sohd, horizontal line, the dom- 
inance of a normal plasma is allowed. F = 15 is assumed through- 
out the jet. 



5 APPLICATION TO THE 3C 279 JET 

Let us next consider the 3C 279 jet. This quasar (^ = 0.538) 
was the first source found to show superluminal motion 
(Cotton et al. 1979; Unwin et al. 1989). At this distance, 
1 mas yr~^ corresponds to /3app = 28.2/i~^. 



5.1 Kinetic luminosity 

In the same manner as in 3C 345 case, we first consider 
dd,coTev/i'o- Since the turnover frequency of the core is re- 
ported to be about 13 GHz (Unwin et al. 1989), we measure 
Sd,core at !/ < 10.7 GHz. At 10.7 GHz, ^d.corc IS reported to 
be 0.79 mas from five-epoch observations by Unwin et al. 
(1989), and 0.61 mas from three-epoch observations by Car- 
rara et al. (1993). With one-epoch VSOP (VLBI Space Ob- 
servatory Programme) observation, 0d,coro = 1.8 x 0.28 mas 
and ^d.corc = 1-8 x 1.95 mas were obtained at 4.8 GHz and 
1.6 GHz, respectively (Piner et al. 2000). As the weighted 
mean over the 10 epochs, we obtain 



— 6.6 X 10 , 

I'D 

which leads to 



(53) 



1 X^TV 

z/0 T\ sin ip 



\ mas / Vo 



3.7 X 10" 



(54) 



Figure 4. The same figure as figure 3, but the different Lorentz 
factor r = 20 is adopted. 



We next consider the spectral index of the core in the 
optical thin frequencies. From mm to IRAS (90-150 GHz and 
270-370 GHz) observations, Grandi et al. (1996) reported 



^1.2. 



(55) 



Adopting this value as the spectral index that reflects the 
energy distribution of the power-law electrons in the core, 
we obtain 



(56) 



/ 1 f^ \ 2{5-2a)/(7-2a 

/ 1 X^TV \ 


= 2.7x ] 


\ i/Q J"! sin ip j 




and 






-laK n' 

7min '0 


-4/(7-2«) 



1 fil V in^l-l r^-0.43„, -0.56 
— 1.81 X 10 A 7min 

If we assume a normal plasma dominance, we obtain 

-1.56 



(57) 



Lkin = 3.0 X 10«ft-^-^^ (:^) 
1 



X r(r-i) 



sin 93 V 1 + 2 



-0.43 



ergs s 



(58) 



5.2 Kinematics near the core 



We now consider V and ^ near to the core to constrain Lkin 
further. The apparent superluminal velocities of the compo- 
nents in the 3C 279 jet do not show the evidence of acceler- 
ation and are roughly constant during the propagation. For 



example, Carrara et al. (1993) presented /3a 



4.51/1" 
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for C3 during 1981 and 1990 over 15 epochs at 5, 11, and 
22 GHz and /3app = 4.23/i"^ for C4 during 1985 and 1990 
over 15 epochs at 11 and 22 GHz. We thus assume that the 
apparent motion of the components represent the fluid bulk 
motion and that both F and <^ are kept constant through- 
out this straight jet. Since /3app does not differ very much 



2.58 ± 0.18 Jy at 22.2 GHz (1984.09) for C3. We modelfit 
the three-frequency data by the function 



between C3 and C4, we adopt /3a 



4.51 as the common 



value and apply the same F and ip for these two components. 
The errors incurred by the difference in /3app are small. 

Let us first consider the lower bound of F. From the ob- 
served /3app, we obtain F > -^//Japp^ -I- 1 = 4.61. As case 1, we 
adopt F = 4.61 and /?app = 4.51. It follows that (/ja = 12.5° 
(or tp = 8.1°) and 5 = 4.6. It is worth comparing this 5 
with that measured in the quiescent state; Unwin derived 
5 = 3.6]|^jg so that the calculated X-ray flux with this S 
agrees with the observed value. Moreover, from the RTXE 
observation in the 0.1-2.0 keV band during 1996.06 and 
1996.11, Lawson and M°Hardy (1998) obtained 2^Jy in the 
quiescent state during 1996.49-1996.54 from 16-epoch obser- 
vations. Applying this X-ray flux density to the spectral in- 
formation {vm, 5m, ce) and 0^ (table 3) for component C4 at 
epochs 1987.40 and 1989.26 and component C3 at 1984.05, 
we obtain 5 > 1.9, > 4.0, and > 3.3, respectively. Therefore, 
we can regard that the value of S — 4.6 are consistent with 
the X-ray observations. In this case, equation (B8[) gives 



Lkin = 6.3x10"^ A' 



-0.57 / 7min \ ' 

V 100 ) 



h ' ergs s (case 1).(59) 



This case will be examined in the next subsection. 

Let us next consider the upper bound of F, which can be 
constrained by the upper bound of 5 for a given /3app ~ 4.51. 
In general, the upper bound of S is difficult to be deduced. In 
the case of 3C 279, however, it is reasonable to suppose that 
a substantial fraction of the X-ray flux in the flaring state 
is produced via synchrotron self Compton (SSC) process. 
Therefore, the Doppler factors ~ 3.9 (Mattox et al. 1993) 
and ~ 5 (Henri et al. 1993) derived for the 1991 flare, and ~ 
6.3 (Wehrle et al. 1998) for the X-ray flare in 1996, give good 
estimates. We thus consider that S does not greatly exceed 
10 for the 3C 279 jet. Imposing S < 10, we obtain F < 6.06 
from /3app ~ 4.61. Combining with the lower bound, we can 
constrain the Lorentz factor in the range 4.6 < F < 6.1. 

As case 2, we adopt a large Lorentz factor of F — 6.0. 
In this case, /3app = 4.51 gives S = 9.8 and ipA = 4.4° (or 
equivalently, ip = 6.8°). There is, in fact, another branch of 
solution that gives smaller Doppler factor of 2.2; however, 
this solution gives too large A''c(SSA)'s (~ lO^cm"^). We 
thus consider only the solution giving S = 9.8 as case 2, 
which results in 

Lkin =4.3x10''^ A-"" fl^^y^'\-^'' ergs s"' (case 2). (60) 

We will consider this case, together with case 1, in the next 
subsection. 



5.3 Composition of individual components 

For the jet component C3, Unwin et al. (1989) presented 
i^m = 6.8 GHz, Sm = 9.4 Jy, a = -1.0, and 8^ = 0.95 mas 
at epoch 1983.1 (table 3). It was also reported in their paper 
that the flux densities were 5.13 ± 0.14 Jy at 5.0 GHz (at 
epoch 1984.25), 4.49 ± 0.17 Jy at 10.7 GHz (1984.10), and 



S„ = Ai[l~- exp(-yl 



22^ 



(61) 



Assuming a = —1.0, we obtain i/m = 6.6 GHz and Sm ~ 5.9 
Jy as the best fit around epoch 1984.10. 

For C4, Carrara et al. (1993) presented fm ~ 11 GHz, 
5*111 ~ 4.3 Jy, a = —0.9, and 8d ~ 0.6 mas from their 1989- 
1990 maps at 5, 11, and 22 GHz. At epoch 1987.4, we use 
the fiux density of 1.43 ± 0.17 Jy at 22 GHz (Carrara et al. 
1993) and that of 3.95 ± 0.20 Jy at 5 GHz (Gabuzda et al. 
1999). We extrapolate the fiux density at 11 GHz from those 
at 1988.17, 1989.26, and 1990.17 presented in Carrara et al. 
(1993) to obtain 3.60±0.20 Jy. Assuming a — —0.9, we can 
model fit the three-frequency data to obtain Um ~ 6.4 GHz 
and 5m = 4.4 Jy. At this epoch (1987.4), component C4 is 
located about 1 mas from the core; therefore, the angular 
size vs. distance relation (Carrara et al. 1993) gives 9d ~ 0.6 
mas. The input parameters for C3 and C4 at these 4 epochs 
are summarized in table 3. 

In case 1 (see § |!|), we obtain iV*(SSA) = 6.4 x 10'' 
and 9.4 x 10^ for C3 at 1983.1 and 1984.10, and 1.2 x 10'' 
and 9.7 x 10^ for C4 at 1987.4 and 1989.26, respectively. 
These large values of A'^(SSA) result in such (unphysically) 
smaU values of iVo (nml)/Ar*(SSA) as 9.6 x 10"^ for C3 at 
1983.1, 1.9 X lO"'^ for C3 at 1984.1, and 1.3 x 10"" for C4 at 
1987.4. Since even a pair-plasma dominance should be ruled 
out when 7V*(nml)/iV*(SSA) < 100 holds (see eq. |||l), we 
consider that the value of 5{= 4.6) is underestimated. We 
thus exclude case 1 from consideration. 

We next examine case 2, in which a larger Doppler 
factor {— 9.8) is adopted. In this case, we obtain reason- 
able values of B and A'^e*(SSA) (table 3). It follows that the 
bestfit values of Af*(nml)/iV*(SSA) becomes 2.8 x 10"^ and 
5.7x 10"^ for C3 at 1983.1 and 1984.10, and 3.2 x 10"^ and 
0.40 for C4 at 1987.4 and 1989.26, respectively. If 5 exceeds 
10 (or equivalently, if F exceeds 6), the density ratio for C4 
at 1989.26 becomes close to unity; therefore, we cannot rule 
out the possibility of a normal-plasma dominance for such 
large Doppler factors. 

On these grounds, we can conclude that the jet com- 
ponents are dominated by a pair plasma with 7min ^ 100, 
provided that 5 < 10 holds. However, if 5 is as small as 5, 
not only a normal-plasma dominance but also a pair-plasma 
dominance is ruled out. Thus, we consider that 5 < 5 is 
unphysical for the 3C 379 jet. 



6 DISCUSSION 

In summary, we derive the proper electron number den- 
sity, A'^^(SSA), of a homogeneous radio-emitting compo- 
nent of which spectral turnover is due to synchrotron self- 
absorption. Comparing A'^(SSA) with the density derived 
from the kinetic luminosity of the jet, we can investigate 
whether the possibility of a normal plasma (e~-p) domi- 
nance can be ruled out or not. Applying this method to the 
3C 345 jet, and using the published spectral data of compo- 
nents C2, C3, C4, C5, and C7, we find that they are domi- 
nated by a pair plasma with 7min ^ 100 at 14 epochs out of 
the total 19 epochs, if we assume F < 15 throughout the jet. 
We also investigate the 3C 279 jet and find that components 
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C3 and C4 are dominated by a pair plasma with 7min <C 100 
at all the four epochs, provided that S < 10. 

It is noteworthy that the kinetic luminosity computed 
from equations (BSl) and (pTI) has, in fact, weak dependence 
on the composition. Substituting Ckin for a pair and a nor- 
mal plasma, we obtain 



Lkin(pair) _ rUc {7- (pair) ) /7min (pair) 
Lkin(nml) mp/27min(nml) 

_ 1 {7-(pair)) 7min(nml) 
9.18 7min (pair) 100 



(62) 



It follows from equation ( [2 7] ) that {7_(pair))/7min(pair) ~ 3 
holds for a = —0.75, for instance. Thus, we obtain compa- 
rable kinetic luminosities irrespectively of the composition 
assumed, if we evaluate them by the method described in 

§1 

Let us examine the assumption of A ~ 1 in equa- 
tions (pd), lp^), and (pd). For a typical value of a = —0.75, 
we obtain ({7_)/7min)-K' = ((5cq/(S) "'"'', where 5cq refers to 
the "equipartition Doppler factor" defined by (Readhead 
1994) 



10-3-6" i?(Q)34 



0.65/1 



-^17-2a 
Z) Uu 



qoz + (go - 1)(-1 + V^qoz + 1) 

Y'^—', (63) 



-35-2a rr 16/, -34 1 l/{13-2ci) 



where F{a) is given in Scott and Readhead (1977); 



Oir 



and 9d are measured in GHz, Jy, and mas, respectively. 
Therefore, we can expect a rough energy equipartition be- 
tween the radiating particles and the magnetic field if Scq is 
close the 5 derived independently. 

Substituting the input parameters presented in table 2, 
we can compute (5cq's for the 3C 345 jet. For component C5 
at 1984.2, we obtain Scq = 736; we thus discard this epoch 
for C5. The computed S^q ranges between 1.1 and 67, with 
the averaged value of 23 and the standard deviation of 9. 
The resultant magnetic field strength, B = 6.7 ± 6.4 mG, 
is reasonable. In the same manner, for the 3C 279 jet, we 
obtain 5 = 9.5 ± 1.5 and 12 ± 8 mG for the four epochs. 

Since (5cq = 23 ± 9 for 3C 345 is comparable with those 
presented in table 2 and since 5 = 9.5 =b 1.5 for 3C 279 is 
comparable with those adopted in § 5.3 (or table 3), we can 
expect a rough energy equipartition in these two blazer jets. 
However, it may be worth noting that the energy density of 
the radiating particles dominates that of the magnetic field 
if & becomes much smaller. Such a case (5 = 4 ~ 12) was 
discussed by Unwin et al. (1992; 1994; 1997), assuming that 
a significant fraction of the X-rays was emitted via SSC from 
the radio-emitting components considered. 

We finally discuss the composition variation along the 
jet. In the application to the two blazers, we assumed con- 
stant Lorentz factors (F) throughout the jet. However, F 
may in general decrease with increasing distance from the 
central engine. If the decrease of F is caused by an entrain- 
ment of the ambient matter (e.g., disk wind) consisting of a 
normal plasma, it follows from equation ( pq ) that 7V^(nml) 
and hence A'^(nml)/A'^(SSA) increases at the place where 
the entrainment occurs. (Note that the kinetic luminosity 
will not decrease due to an entrainment in a stationary jet.) 
Therefore, the change of composition from a pair plasma 



into a normal plasma along the jet may be found by the 
method described in this paper. 

To study further the change of the composition, we need 
detailed spectral information (i^m, ^m, a) of individual com- 
ponents along the jet. To constrain the spectral turnover ac- 
curately, we must modelfit and decompose the VLBI images 
from high to low frequencies; in another word, spatial reso- 
lution at low frequencies is crucial. Therefore, simultaneous 
observations with ground telescopes at higher frequencies 
and with space+ground telescopes at lower frequencies are 
preferable in this study. 
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APPENDIX A: MAGNETIC FIELD STRENGTH 
DEDUCED FROM SYNCHROTRON 
SELF-ABSORPTION THORY 

The magnetic field strength, B, can be computed from the 
synchrotron spectrum, which is characterized by the peak 
frequency, fm, the peak fiux density, S^, and the spec- 
tral index, a, in the optically thin regime. Marscher (1983) 
considered a uniform spherical synchrotron source and re- 
lated B with i/ui, Sm, and a in the optically thin cases 
(0 > a > —1.0). Later, Cohen (1985) considered a uniform 
slab of plasma as the synchrotron source and derived smaller 
values of B for the same set of {i/jSmja). In this appendix, 
we express B in terms of (!/,Sm,a) for a uniform sphere for 
arbitrary optical thickness (i.e., a < 0). We assume that the 
magnetic field is uniform and that the distribution of radiat- 
ing particles are uniform in the configuration space and are 
isotropic and power-law (eq.fel]) in the momentum space. 

For a uniform B and A^'^ , the transfer equation gives 
the specific intensity 



I,,* — Au* ' [1 — exp(— a^*a;o*)] . 



where 



Aia) 



(I)" 



e a{a) / e 
c C{a) \27rmcC 



-3/2 



(Al) 



(A2) 



and xo* gives the physical thickness of the emitting region 
along the li ne o f sight. At the distance 6 away from the cloud 
center (fig. Al), the fractional thickness becomes 



2R' 



R^ 



402 



9r 



(A3) 



which is Lorentz invariant. Thus, the specific intensity at 6, 
or equivalently at angle 6 away from the center, becomes 
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Figure Al. Flux emitted from a spherical cloud of radius R* 



lAO) 



l + z 



1/2 



Av 



'/2 



exp 



402 

--40),/!-^ 



,(A4) 



where r^(0) = a^* ■ 2R* is the optical depth for 6 = Q (or 
6 = 0); we utilized Ivjv''' — Iv* jv*'^ and equation (W). Even 
if special relativistic effects are important, the observer find 
the shape to be circular. We thus integrate /^ (ff) cos Q over 
the emitting solid angles 2-k sin QdQ in < < 0d/2 to obtain 
the flux density 



S^ = 2-K 



= "L (hy 



ed/2 



/,.(6')cos6'sin6'de 



4 V rad / V 1 + z 



1/2 



Av 



5/2 



(A5) 



where x = 49^ /dd'^ and 61 < 1 is used. 

Putting the derivative dSv/du to be 0, we obtain the 
equation that relates r^ (0) and a at the turnover frequency. 



-T„(0)vT^^ 



dx 



(l - ^a) r r.(0)yT^e--(o'-^dx. 



(A6) 



We denote the solution r^(0) aX v = Vr^ as rm(0), which is 
presented as a function of a in table 1. 



It is worth comparing equation (A5) with that for a 
uniform slab of plasma with physical thickness R* . If the 
slab extends over a solid angle 7r(0d/rad)^/4, the flux density 
becomes 



^"^ 4 (rad) (l + J 



1/2 



Av 



.V2/l_e-r„(0)^ 



(A7) 



The optical depth Tv{Q) — at ■ 2R* does not depend on h 
(or 6) for the sl ab g eometry. Therefore, comparing equa- 
tions (A5) and (A_7,) we can deflne the effective optical 
depth, (r^), for a uniform sphere of plasma as 



-<r„) 



-r^(0)VT^ 



dx- 



(A8) 



We denote {t^) at i^ = j/m as (rm), which is presented in 
table 1. Note that (rm) is less than rm(0) due to the geo- 
metrical factor -yi — X- 

Using rm(0), we can evaluate the peak flux density, Sm, 
at V = Vm and inversely solve equation (A.5) for B. We thus 
obtain equation (Ud) with 



6(q) 



X 10' 



(1)" 



a{a) 



C{a) 
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„(0)^^^^ 



dx 



(A9) 



The values are tabulated in table 1; they are, in fact, close 
to the values obtained by Cohen (1985), who presented 6(a) 
for a slab geometry with —a < 1.25. This is because the 
averaged optical depth, (r^) at i^ = j/m, for a spherical ge- 
ometry becomes comparable with that for a slab geometry 
(Scott & Readhead 1977). 

If we took the optically thin limit, rm(0) <C 1, we could 
expand the integrant in equation (A5) as 



-T„,(0)VT^O? 



.(0)v/l^, 



(AlO) 



to obtain fe(— 0.5) — 3.33 and 6(— 1.0) = 3.83, for example. 
This case (rm <C 1) has been considered in many papers for 
smaller —a (< 1.0). 
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Table 2: Input parameters and the results for the 3C 345 jet when F = 15 is assumed. 



epoch 


p sin (^t 
mas 


/3app/lt 


GHz 


Si 
Jy 


at 


mas 


5i 


deg 


mG (68 % errors) 


cm"^ 


7V;(nml)/7V;(SSA)t 
(68 % errors) 














component 


C2 










1982.00^ 


4.9 


12.9 


1.5 


2.0 


-0.60 


2.15 ±0.18 


22.5 


2.19 


13.61 ( 9.42, 18.36) 


0.98E-3 


0.85E±0 ( 0.44E±0, 0.15E±1) 














component 


C3 










1981.50'' 
1982.00^ 


2.0 
2.2 


8.7 
9.2 


2.7 ±0.3 
2.6 


2.4 ±0.2 
2.1 


-0.80 
-0.70 


1.00 
0.97 ±0.07 


27.2 
26.8 


1.22 
1.32 


9.74 ( 5.23, 16.54) 
9.46 ( 6.85, 12.18) 


0.47E-1 
0.24E-1 


0.82E-1 ( 0.29E-1, 0.21E±0) 
0.17E±0 ( 0.95E-1, 0.29E±0) 














component 


C4 











1982.00^ 0.4 6.2 14.6 7.6 -0.30 0.29 ± 0.02 28.6 0.82 26.90(19.94,33.87) 0.34E-1 0.13E±1 

1983.40'' 0.7 6.0 11.6 ±0.4 7.5 ± 0.6 -0.50 0.37 28.7 0.81 29.05(22.55,35.74) 0.64E-1 0.44E±0 

1985.80'' 1.0 7.9 9.4 ± 0.6 1.9 ±0.1 -0.20 0.53 27.7 1.09 386.49(273.61,521.36) 0.36E-4 0.38E±3 

1988.20^ 2.0 14.2 6.5 ± 1.4 0.8 ±0.1 -0.60 0.69 19.7 2.75 1206.26(347.89,3172.78) 0.39E-4 0.21E±3 

1990. 70<= 4.3 14.0 2.0 ± 0.5 2.0 ± 0.5 -0.75 ±0.15 0.85 20.3 2.63 1.24(0.28,4.20) 0.61E±1 0.87E-3 

1992.50'' 4.3 13.4 2.7 ±1.6 1.6 ±0.1 -0.70 0.97 21.7 2.36 15.91(0.87,209.80) 0.24E-1 0.17E±0 



0.84E±0, 0.20E±1) 
0.29E±0, 0.65E±0) 
0.24E±3, 0.54E±3) 
0.32E±2, 0.95E±3) 
0.46E-4, 0.13E-1) 
0.50E-2, 0.19E±2) 



component C5 



1984.20" O 5i) 2.5 ±0.4 4.9 ±1.1 -0.80 0.27 29.1 0.66 0.01 ( 0.00, 0.02) 0.59E+6 0.89E-7 

1985.80'' 0.6 6.4 6.0 ± 0.3 3.3 ± 0.3 -0.60 0.37 28.5 0.86 5.68(4.11,7.49) 0.63E±0 0.45E-1 

1987.30'' 0.9 7.5 7.0 ± 0.6 5.1 ± 0.2 -0.80 0.47 28.0 1.02 12.68(8.00,18.75) 0.99E±0 0.18E-1 

1988.20'' 1.1 8.0 7.5 ±0.4 5.6 ± 0.5 -0.70 0.53 27.6 1.11 24.42(17.69,32.94) 0.12E±0 0.11E±0 

1989.20'' 1.4 8.6 2.5 ± 1.4 4.7 ±0.6 -0.90 0.59 27.2 1.21 0.20(0.01,2.14) 0.14E±4 0.76E-5 

1990.30'' 1.7 9.3 2.5 ± 1.2 3.6 ±0.1 -0.90 0.66 26.7 1.33 0.52(0.04,4.07) 0.15E±3 0.60E-4 

1990. 55^ 1.8 9.5 2.7 ±0.5 3.2 ± 0.5 -0.75 0.80 26.5 1.37 2.23(0.77,5.47) 0.13E±1 0.46E-2 

1990. 70<= 1.7 11.3 2.7 ±1.5 3.2 ± 0.9 -0.75 ±0.15 0.80 ± 0.20 24.8 1.74 2.08(0.09,30.81) 0.22E±1 0.27E-2 

1992.50'' 2.3 11.2 2.2 ± 0.8 1.5 ± 0.3 -0.90 0.81 24.9 1.72 3.37(0.40,17.30) 0.14E±1 0.43E-2 



0.16E-7 
0.26E-1 
0.81E-2 
0.60E-1 
0.21E-6 
0.17E-5 
0.80E-3 
0.50E-4, 
0.16E-3, 



0.48E-6) 
0.74E-1) 
0.34E-1) 
0.19E±0) 
0.12E-2) 
0.38E-2) 
0.22E-1) 
0.14E±1) 
0.11E±0) 



component C7 



1992.05^ (U 2^8 12.8 ± 0.5 4.6 ± 0.5 -0.75 0.20 ± 0.04 29.7 0.36 11.29(4.39,24.14) 0.43E±1 

1992.67'' 0.2 6.0 12.5 ±1.0 7.0 ± 0.5 -0.75 0.35 ± 0.02 28.7 0.80 39.25(23.65,60.79) 0.33E±0 

1993.19'' 0.4 10.5 11.6 ±0.5 5.1 ± 0.5 -0.75 0.41 ± 0.02 25.7 1.55 85.79(59.35,118.42) 0.96E-1 

1993.55'' 0.5 14.5 11.0 ±1.5 3.1 ± 1.0 -0.75 0.38 ± 0.02 20.9 2.56 106.85(41.57,294.30) 0.14E±0 



0.22E-1 
0.94E-1 

0.24E±0 



0.35E-2, 0.10E±0) 
0.38E-1, 0.21E±0) 
0.12E±0, 0.45E±0) 



0.18E±0 ( 0.30E-1, 0.14E±1) 



For the input parameters, 1 a errors are presented if they are given in the litearture. 

For the output parameters, 1 a errors are presented in the parentheses for some parameters. 

Input parameters are cited from Zensus et al. (1995). 

Input parameters are cited from Lobanov and Zensus (1999). 

Input parameters are cited from Unwin et al. (1994). 

Input parameters are cited from Unwin et al. (1997). 



Table 3: Input parameters and the results for the 3C 279 jet in case 2. 



epoch 


psini^ 
mas 


Papp*^ 


GHz 


Jy 


a 


mas 


5 


<y3A 

deg 


B 
mG 


^e(SSA) 
cm~^ 


N*{ri^\)lN*{S?>K) 














component 


C3 










1983.10^ 
1984.10'' 


1.4 
1.4 


4.5 
4.5 


6.8 
6.5 


9.4 

8.9 


-1.00 
-1.00 


0.95 
0.80 


9.8 
9.8 


4.45 
4.45 


17.67 

8.24 


0.54E-I-2 
0.38E-F3 


0.28E-2 
0.57E-3 














component 


C4 










1987.40" 
1989. 26<= 


1.0 
1.2 


4.2 
4.2 


6.4 
11.0 


4.4 
4.3 


-0.90 
-0.90 


0.60 
0.60 


9.8 
9.8 


4.45 

4.45 


9.80 
157.89 


0.12E-I-3 
0.94E-F0 


0.32E-2 
0.40E-I-0 



* Input parameters are cited from Unwin et al. (1989). 
'' Input parameters are computed in this paper (see text) 
'^ Input parameters are cited from Carrara et al. (1993). 



